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Introduction 

Radiation transport models are required for the 
evaluation of risk to passengers and crews on high- 
altitude airplanes and to astronauts in low-Earth or- 
bit or deep space and for understanding degradation 
properties of materials, including microelectronics, 
on orbiting satellites. Light ions such as 2 H or 4 He 
are present as primary and secondary radiations in 
solar and galactic cosmic rays. Light ion beams also 
are useful in radiotherapy treatment of cancer pa- 
tients. (See refs. 1 and 2.) In nuclear reactions the 
mass dependence of secondary radiations results in 
broader distributions of secondary energies for lighter 
mass particles. Several mechanisms (direct knock- 
out, evaporation, coalescence, etc.) will produce light 
products of nuclear collisions and lead to a wider dis- 
tribution of secondary energies. The broad spectrum 
of light ion secondaries leads to some different ap- 
proaches to solution of the Boltzmann equation for 
radiation transport compared with heavy ion trans- 
port. (See ref. 3.) The heavy fragments in nuclear 
reactions are usually spectators in that they receive 
only small momentum and energy transfer. 

In this report, we discuss the modifications of 
the monoenergetic version of the BRYNTRN trans- 
port code (refs. 4-6) for the transport of light ion 
beams. Validation of radiation transport methods 
must include monoenergetic beams and a continu- 
ous spectrum of primary radiations. The BRYNTRN 
transport code solves the Boltzmann equation in 
the straight-ahead approximation using the method 
of characteristics developed in references 3-6. Re- 
cently we developed an energy-dependent model of 
light ion nuclear interaction cross sections. (See 
refs. 7 and 8.) The previous monoenergetic ver- 
sion of the BRYNTRN code considered only the cou- 
pled proton-neutron transport problem. Here we dis- 
cuss the extension of the BRYNTRN code to include 
the database for transport of the light ion beams 
2 H (d), 3 H(i), 3 He(h), and 4 He(a). The BRYNTRN 
code previously calculated target fragmentation ef- 
fects through local energy deposition. We also de- 
scribe an equilibrium target fragment field that was 
added to consider the effects of fragments on the lin- 
ear energy transfer (LET) spectrum. 

Symbols 

A mass number 

A? mass number of target 

c speed of light, fm/sec 

D dose, cGy 


d 2 H 

E particle energy, MeV/amu 

E' = £ (r') 

Eq beam kinetic energy, MeV /amu 

Eqj average energy of target fragment j 

F defined in equation (15) 

F defined in equation (14) 

fjk(E, E r ) differential energy cross section for 
redistribution of particle type and 
energy 

h 3 He 

m nucleon mass, MeV/c 2 

R nuclear matter radii, fm 

r residual range, g/cm 2 

S(E) proton-stopping power, MeV/(g/cm 2 ) 

t 3 H 


x 

Z 
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a 
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fijk 

s(r) 
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a j(E) 
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Subscripts: 


depth of material, g/cm 2 
charge 

distance, g/cm 2 
alpha particle 
Dirac delta function 
Kronecker delta 

energy associated with r, MeV /amu 

defined in equation (7) 

range-scaling parameter 

cross section, mb 

total cross section, mb 

linear energy transfer, MeV/cm 

differential flux spectrum of type j 
ions 

scaled flux 


EL elastic recoils 

j , k ions of type j and k 

n neutron 

T target 



Boltzmann Equations and Solutions 

The propagation of high-energy nuclei and their secondaries through bulk matter is described by the 
Boltzmann equation, which in the straight-ahead approximation (refs. 5 and 6) is of the form 


9 v j’§E S ^ + a A B ) 


dx 


<i>j(x,E) = Ef f jk {E,E')ct> k (x,E')dE' 


k Jo 


( 1 ) 


where Vj denotes the range-scaling parameter equal to Zj/Aj, Z is the charge, and A is the mass number. In 
equation (1), 5(E) is the proton-stopping power, <j(E) is the total cross section, ^(x, E) is the differential flux 

spectrum of type j ions, and f jk (E, E') is a differential energy cross section for redistribution of particle type 
and energy. 


Utilizing the definitions 


and 

we can rewrite equation (1) as 
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r ~Jo 5(E0 
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dx 3 dr 

which is solved by Wilson et al. (refs. 5 and 6) as 

i/jj(x, r) = exp[— ^-(r, x)]ipj (0, r + v jX ) 

+ [ n / ex p[~c j( r , z )\fjk( r + l/ jZ,r , )'(p k (x — z,r')dr , dz 
k Jr 

where the exponential is the integrating factor given by 

Cj (u 0 ~ J <Xj (r + dt' 


(2) 

( 3 ) 

( 4 ) 

(5) 
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Next we present the basis of numerical procedures for propagation of the solution at ip Ax, r) to tyAx + hr ). 
Choosing h to be small such that 

<jj(r r )h < 1 

then from perturbation theory (ref. 4) we derive 


( 8 ) 

( 9 ) 


^(x + h-z, r') « exp [~( k ( r',h - z)\ ^ k [x, r' + v k {h- z)\ 

For monoenergetic particle beams, the initial beam of type J particles of energy E 0 is 

^■(0, r) = 8jj8(rQ - r) (10) 

where rg — R(Eq). The separation of the singular terms in the solution is written by the replacement 

ipj(x,r) -^'ip j o(x,r) + ijj j (x,r) 
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( 11 ) 



Based on equations (8)— (11), the solution (refs. 3 and 7) at x + h is 
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for charged-particle propagation, where r' Qj = r 0 - Uj[x + (h/ 2)], and 
rp n {x + h, r) « exp[-£ n (r, /i)]^ n (x, r) 


+*E( exp<J — 

k 


Cn( r ,^) + a(r^ 




+ Jr dr' exp | Cn (V 2 ) + & }/n* (g r 0^ (z, r f + i/*. ^ 


for neutron propagation. In equation (12) F is related to the cumulative spectrum F as 


(13) 


fh _ 

Fij (r, Vjh, r') — / /y (r + z'jz, r') dz = F^ (r + i/jh, /) - F t j (r, r') 

J 0 


(14) 


with 



fij {E, E r ) dE 


(15) 


where e(r) is the energy associated with the residual range r and E' — e(r'). 

The numerical representation of particle ranges and stopping powers in the BRYNTRN code has been 
described in reference 6, where the database for proton and neutron nuclear interactions is also described. The 
database for light ion interactions is described in references 8 and 9. We note that the projectile-like fragments 
are represented by a Gaussian spectrum with energy-dependent width and downshift from the beam velocity. 
We expect that a sum of several Gaussian terms would more accurately represent these spectra. (See ref. 10.) 
The quasi-elastic events for light ions are represented by an approximation to the first collision term of the 
inelastic scattering series. (See ref. 11.) 


Secondary Flux From Light Ion Beams 

In figures 1 and 2, we show the secondary charged particle spectra at several depths in water for 
monoenergetic He beams at 200 and 1000 MeV/amu, respectively. Two sources of secondary particles are 
clearly seen in the figures that correspond to a high-energy flux from projectile-like fragments and a low-energy 
flux from target-like fragments. The projectile-like fragments are produced at energies above that of the beam 
due to the internal Fermi motion of the nuclear-bound state. For the beam at 200 MeV/amu and for great 
depths, the projectile- and target-like secondaries overlap as the beam energy degrades and the secondaries 
slow from atomic interactions. A similar comparison is shown in figures 3 and 4 for aluminum shielding. In 
figures 1-4, the greater range of 3 H compared with 3 He is apparent. 

In figures 5 and 6, we show the secondary proton spectra for monoenergetic 2 H beams at several depths 
in aluminum for 200 and 1000 MeV/amu, respectively. The pickup reactions for 2 H projectiles have not been 
included at this time, and no projectile-like secondaries with £ > 1 are included. Before discussing depth- 
dose curves and the LET spectrum for light ion beams, we next describe an equilibrium field of heavy target 
fragments of mass number Ap >4. 
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Equilibrium Field of Target Fragments 

The nuclear secondaries produced from the target atoms with Ap > 4 have average energies of only a few 
MeV or less and, therefore, are of extremely short range. Previously, their effects had been included only as 
locally deposited energy. (See ref. 6.) Consideration of the target fragmentation effects in the LET spectrum 
and of the radiation response quantities other than energy deposition requires an equilibrium field of target frag- 
ments. The differential flux for these ions can be solved in closed form (ref. 6) as 

1 f°° 

M x < E '< E k)= g-^}2J E fkj(E',E k )4> k (x,E k )dE' (16) 

where <p k (x, E k ) is a projectile-like ion k with energy E k at x, which produces a target fragment j with energy E. 
The production spectrum from fragmentation is assumed of the form 


where p T is the density of the target material, aj is the microscopic fragmentation cross section, and E 0 j is 
the average target fragment energy. The parameters Oj and E 0j used in the BRYNTRN code are described in 
reference 6. Using equation (17) in equation (16), we find 


E \ 2 [ E ( —E V 

2B 0 J ^V2£o/ XP (2^Jj (18) 

where Erfc is the complementary error function. For composite materials, a summation over the target 
constituents is implied in equation (17). The energy spectrum of elastic recoils (ref. 6) is approximated by 


(j>j{x,E]E k ) = 


j{ E k)PT 

Sj (E) 


4>k(E k ,x)\Erfc\ 


fkj( E > E k) = &ELPT 


Cel exp(— C el E) 

1 - exp[-C EL (l - D EL )E k ] 


(19) 


where a EL is the elastic scattering cross section and C EL = A T mc 2 [B + (a 3 /3)], where A T is the mass number 
of the target, m is the nucleon mass, B is an average nucleon-nucleon slope parameter, and a is the target 
radius parameter. The energy transfer to the recoil is restricted by 

0 < E < (1 - D Eh )E k (20) 

with 


The flux of recoils is then 


d EL = 


{A t - 1) 2 
{At + l ) 2 


( 21 ) 


(f)j{x,E;E k ) 


PT a EL { exp ( — C ee E) - exp[-C EL (l - D EL )E k }} , 

S j( E ) 1 - exp[-C EL (f- D EL )E k ] VkW, k) 


The equilibrium field of target fragments is found as 


( 22 ) 


0 , E ) — J (f)j {. X , E] E k ) dE k (23) 

For calculations, the production cross sections from light ions are assumed to scale to the nucleon-induced ones 
by the factor A° k A . (See ref. 3.) An energy grid of 100 points with a maximum energy of 30 MeV/amu is used 
to obtain convergence of the equilibrium target fragment fields for water. Interface effects at the boundary of 
two media could be described in a similar manner as followed here and as described in reference 12. 
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Depth-Dose Curves for Monoenergetic 
Beams 

In figures 7-9, we show calculations of absorbed 
dose versus depth in water shields for monoenergetic 
light ion beams at several energies. These figures pro- 
vide a breakdown of the contributions from primaries 
(dotted line), secondary protons (dash-dot line), and 
all other secondaries (dashed line). Included in the 
secondary contribution are the contributions from 
heavier target fragments (Ap > 4) as described ear- 
lier. The increase in secondary contributions to the 
absorbed dose is apparent when the proton beam is 
compared with the heavier beams where projectile 
fragmentation provides additional particles. The 2 H 
range is about twice that of the proton; therefore, nu- 
clear reactions can continue, which is also apparent 
in figures 7-9. The absorbed dose does not effectively 
represent harmful biological components of a radia- 
tion field, so we next consider the LET spectrum, 
which is more useful. 

LET Spectrum for Monoenergetic 
Beams 

The LET spectrum is closely related to the re- 
sponse of microdosimeters and provides a rough in- 
dication of biologically harmful components of a 
radiation field. Wilson and Badavi (ref. 13) have dis- 
cussed the numerical treatment of the singularities in 
the transformation of energy spectra to LET spectra. 
We use this method to consider the LET spectrum 
in water for monoenergetic light ion beams and all 
nuclear secondaries from target fragmentation. 

In figure 10(a), the LET spectrum from a 
200-MeV proton beam is shown at various depths 
in water. A similar comparison is provided in fig- 
ures 10(b) and 10(c) for a 200 MeV/amu, 2 H and 
4 He beam, respectively. The inclusion of the heavy 
component of the target fragments (Ap > 4) pro- 
vides particles above 200 keV / pm with significant 
numbers at the greater depths. The spikes in fig- 
ures 10 and 11 occur both at a LET value corre- 
sponding to the peak in a Bragg ionization curve 
and near the LET value corresponding to minimal 
ionization for each charge group. The biological ef- 
fectiveness of these ions is not a unique function of 
LET because biological effectiveness also depends on 
ion velocity due to thindown effects where the range 
of secondary electrons becomes small compared with 
the size of relevant damage sites inside cells. (See 
ref. 14.) In figures ll(a)-ll(c), calculations of the 
LET flux spectrum at 1000 MeV/amu are shown. 
The high LET components increase compared with 
the 200 MeV/amu beams in figures 10(a)-10(c) due 
to more predominant nuclear fragmentation effects. 


A similar comparison for the integral LET spectra is 
shown in figures 12(a)-12(c) and 13(a)-13(c). 

Conclusion 

The monoenergetic version of the BRYNTRN 
transport computer code has been extended to the 
transport of light ion beams ( 2 H, 3 H, 3 He, and 4 He). 
Illustrative calculations of energy spectrum, depth- 
dose curves, and LET spectrum were presented and 
discussed. The extensions described here will allow 
transport methods for space radiations to be vali- 
dated through comparisons with laboratory experi- 
ments that involve light ion beams. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
August 16, 1994 
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Figure 1. Calculations of light ion flux spectrum in water shields from primary 4 He beam at 200 MeV/amu. 
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(a) Alpha emissions. 


10~ 3 
7^ 

3 

E 
ro 

a> 10 -4 

I 

C\J 

E 
(5 

•I 10- 5 

Q_ 

v— 

0) 

CL 
X 

LL 

10~6 
10 

Energy, MeV/amu 

(b) Helion emissions. 

Figure 2. Calculations of light ion flux spectrum in water shields from primary 4 He beam at 1000 MeV/amu. 
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Figure 3. Calculations of light ion flux spectrum in aluminum shields for primary 4 He beam at 200 MeV/amu. 
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Figure 4. Calculations of light ion flux spectrum in aluminum shields for primary 4 He beam at 1000 MeV/amu. 1 
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Figure 5. Calculations of proton flux spectrum in water shields from primary 2 H beam at 200 MeV/amu. 
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Figure 6. Calculations of proton flux spectrum in water shields from primary 2 H beam at 1000 MeV/amu. 
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Proton beam. 
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Figure 8. Depth-dose curves for beams at 200 iVleV j amu in water. 
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Figure 13. Integral LET flux for beams at 1000 MeV/amu and several depths in water. 


0(>LET) per primary, cm -2 



<5(>LET) per primary, cm -2 


T) 

o 

C“t“ 

o 

3 

cr 

CD 

P 


m 

H 

7T 

<D 

< 

3 






REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


«th«int P anH m^nlafnin^ ri? rt ^ ^ f ln f orrnatl °" ls est| m ated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and review.ng the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports 1215 Jefferson 
^ ton 1 VA 22202 ~ 4302 ' and to the 0ffice of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 


1 AGENCY USE ONLYpeave blank) 2. REPORT DATE 

November 1994 


4. TITLE AND SUBTITLE 

Extension of the BRYNTRN Code to Monoenergetic Light Ion 
Beams 


3. REPORT TYPE AND DATES COVERED 

Technical Paper 


5. FUNDING NUMBERS 

WU 199-45“ 16- 11 


6. AUTHOR(S) 

Francis A. Cucinotta, John W. Wilson, and Francis F. Badavi 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

NASA Langley Research Center 
Hampton, VA 23681-0001 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 

L- 17394 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA TP-3472 


11. SUPPLEMENTARY NOTES 

Cucinotta and Wilson: Langley Research Center, Hampton, VA; Badavi: Christopher Newport University 
Newport News, VA. J 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified-Unlimited 


12b. DISTRIBUTION CODE 


Subject Category 93 


13. ABSTRACT (Maximum 200 words) 

A monoenergetic version of the BRYNTRN transport code is extended to beam transport of light ions ( 2 H, 
H, He, and He) in shielding materials (thick targets). The redistribution of energy in nuclear reactions is 
included in transport solutions that use nuclear fragmentation models. We also consider an equilibrium target- 
fragment spectrum for nuclei with mass number greater than four to include target fragmentation effects in 
the linear energy transfer (LET) spectrum. Illustrative results for water and aluminum shielding, including 
energy and LET spectra, are discussed for high-energy beams of 2 H and 4 He. 


14. SUBJECT TERMS 

Radiation transport; Cosmic rays; Radiotherapy 


15. NUMBER OF PAGES 

34 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


NSN 7540-01-280-5500 


16. PRICE CODE 

A03 


18. SECURITY CLASSIFICATION 19. SECURITY CLASSIFICATION 20. LIMITATION 
OF THIS PAGE OF ABSTRACT OF ABSTRACT 

Unclassified Unclassified 


Standard Form 298(Rev. 2-89 


Prescribed by ANSI Std. 239-18 
298-102 














